home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zlaic1.z / zlaic1
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZLLLLAAAAIIIICCCC1111((((3333FFFF))))                                                          ZZZZLLLLAAAAIIIICCCC1111((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZLAIC1 - applie one step of incremental condition estimation in its
  10.      simplest version
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZLAIC1( JOB, J, X, SEST, W, GAMMA, SESTPR, S, C )
  14.  
  15.          INTEGER        J, JOB
  16.  
  17.          DOUBLE         PRECISION SEST, SESTPR
  18.  
  19.          COMPLEX*16     C, GAMMA, S
  20.  
  21.          COMPLEX*16     W( J ), X( J )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      ZLAIC1 applies one step of incremental condition estimation in its
  25.      simplest version:
  26.  
  27.      Let x, twonorm(x) = 1, be an approximate singular vector of an j-by-j
  28.      lower triangular matrix L, such that
  29.               twonorm(L*x) = sest
  30.      Then ZLAIC1 computes sestpr, s, c such that
  31.      the vector
  32.                      [ s*x ]
  33.               xhat = [  c  ]
  34.      is an approximate singular vector of
  35.                      [ L     0  ]
  36.               Lhat = [ w' gamma ]
  37.      in the sense that
  38.               twonorm(Lhat*xhat) = sestpr.
  39.  
  40.      Depending on JOB, an estimate for the largest or smallest singular value
  41.      is computed.
  42.  
  43.      Note that [s c]' and sestpr**2 is an eigenpair of the system
  44.  
  45.          diag(sest*sest, 0) + [alpha  gamma] * [ conjg(alpha) ]
  46.                                                [ conjg(gamma) ]
  47.  
  48.      where  alpha =  conjg(x)'*w.
  49.  
  50.  
  51. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  52.      JOB     (input) INTEGER
  53.              = 1: an estimate for the largest singular value is computed.
  54.              = 2: an estimate for the smallest singular value is computed.
  55.  
  56.      J       (input) INTEGER
  57.              Length of X and W
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZLLLLAAAAIIIICCCC1111((((3333FFFF))))                                                          ZZZZLLLLAAAAIIIICCCC1111((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      X       (input) COMPLEX*16 array, dimension (J)
  75.              The j-vector x.
  76.  
  77.      SEST    (input) DOUBLE PRECISION
  78.              Estimated singular value of j by j matrix L
  79.  
  80.      W       (input) COMPLEX*16 array, dimension (J)
  81.              The j-vector w.
  82.  
  83.      GAMMA   (input) COMPLEX*16
  84.              The diagonal element gamma.
  85.  
  86.      SEDTPR  (output) DOUBLE PRECISION
  87.              Estimated singular value of (j+1) by (j+1) matrix Lhat.
  88.  
  89.      S       (output) COMPLEX*16
  90.              Sine needed in forming xhat.
  91.  
  92.      C       (output) COMPLEX*16
  93.              Cosine needed in forming xhat.
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.